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ABSTRACT 

We present continuous, high-precision photometric monitoring data with 1 
minute cadence of the dM3e flare star AD Leo with the MOST satellite. We 
observed 19 flares in 5.8 days, and find a flare frequency distribution that is 
similar to previous studies. The light curve reveals a sinusoidal modulation with 
period of 2.231qJ7 days that we attribute to the rotation of a stellar spot rotating 
into and out of view. We see no correlation between the occurrence of flares and 



1 Based on observations obtained with the Apache Point Observatory 3.5-meter telescope, which is owned 
and operated by the Astrophysical Research Consortium. 
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rotational phase, indicating that there may be many spots distributed at different 
longitudes, or possibly that the modulation is caused by varying surface coverage 
of a large polar spot that is viewed nearly pole-on. The data show no correlation 
between flare energy and the time since the previous flare. We use these results 
to reject a simple model in which all magnetic energy is stored in one active 
region and released only during flares. 



Introduction 



Flares are explosive events caused by magnetic reconnection in stellar atmospheres. On 
M dwarfs, where flare signat ures can be spectacular, they emit at a range o f wavelengths, 
from radio JOsten et alJbood) and I R fl^midtll20nypavenportl[2^1^ to UV (Uudard et aL 



2000; Sanz-Forcada fc Micela 2002; Welsh et al. 



2007h and X-ray (fasten et alJboiol 



tailed studies of flares, especially m ulti-wavelength observations (e.g., lHawley et al. 



De- 



2003 



Osten et al.ll2005l ; ISmith et al.l 120051 ) . allow for interpretation of th e physical processes that 
occur during flares (for a thorough review, see lBenz fc Gudelll2010l ). 



The emission from flares can also be exploited as an indicator of the strength and 
distribution of the magnetically active regions on the stellar surface. The frequency of flares 
of a given energy, or the flare frequency distribution (FFD), of a star provides one such 
metric by tracking how much magnetic energy is released through flares. Photometric flare 
monitoring efforts (e.g.. lGershberglll972l ; lLacy et al.lll976l ). typically undertaken in Johnson 
U, have found that more energetic flares occur less frequently, according to 



log v = a + (3 log E 



where v is the frequency of flares of energy E or greater, and j3 is the power law exponent, 
typically < -1. 

Studies of the FFD of a number of stars have found that later- typ e M dwarfs flare more 



frequently but wi t h less energy than their earlier- type counterparts (IKowalski et al.l 12010 



Hilton et al. 


2010; 


Hilton 


2012) 



an indicator of magnetic activity, flare more frequently and with more energy than inactive 
stars. 



The FFD for the subi 



several times. lLacy et al. 



e ct of t his study, the bright dM3e star AD Leo, has been measured 
(119761 ) found (3 = —0.82 ± 0.27, although this measurement was 



based on only 9 flares in 21. 5 hours of monitoring. With extensive monitoring over several 



years, 



Pettersen et al.l (11984 ) found a value of j3 = — 0.62±0.09, and did not see differences 
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above their detection limits in the FFD from year to year, thus arguing for no cyclic variations 
caused by a magnetic cycle. They also found that flares were distributed randomly in 
time, inferring that flares were occurring either from one spot always in view (nearly pole- 
on) or from ma ny spots. The form er scenario is further supported by spectropolarimetric 



observations of iMorin et al.l ( 120081 ). who report that the dominant feature of the surface 
magnetic field is likely a large polar spot, vi ewed nearly pole on , with radial field of maximum 
magnetic flux B = 1.3kG. We note that the lMorin et al.l ( 120081 ) result is based only on Stokes 



kG ( 


Johns-Kru 


11 & Valenti 


2000 


), \Bf\ 


(Shulvak et al.l 


2010 


)• 



2.9 kG ( IReiners fc Basril 120071 ). and |5| = 2 



3.3 
2.5 



A new generation of space missions allows for previously inaccessible continuous, short 
cadence, high precision photometric monitoring of flares. These satellites have sufficient 



precision to measure rotational modula tion caused by star spots (e.g., iRucinski et al.l 12004 ; 
Basri et al.l 1201 ll ; lHarrison et al.l 120121 ) and to detect small flares. Practically continuous 
monitoring ensures that nearly all flares are observed. For the first time, we can use the 
relative timing of flare occurrence to test simple models of the distribution of the active 
regions on M dwarfs. 

We present the results of 8 days of continuous photometric monitoring observations of 
AD Leo with the Microvariability and Oscillations of STars (MOST) satellite. The details 
of the observations and photometry are provided in Section 121 along with our methods for 
processing the light curve and for identifying flares. Our measurements of the FFD and an 
analysis of the timing of the observed flares are described in Section [3j where we compare 
these results to a simple, single active region model. We conclude in Section [U 



2. MOST Photometry 

The MOST mission is a Canadian microsatellite tasked with the detection and charac- 
terization of acoustic oscillations in Sun-like stars, metal-poor subdwarfs, and magnetic stars 
in order to seismically probe their structures and ages. MOST travels in a sun-synchronous, 
polar orbit with a period of 101.4 minutes, allowing it to observe a single targe t within the 



Contin uous Viewing Zone for up to ~60 days at a time without interruption ( iRowe et al. 



2006aJ). The sensitivity and speed of MOST provides micromagnitude preci sion for stars 



with 0.4 < V < 5.5 with individual exposure times of less than 1 minute (I Walker et al. 



20031 ). The less precise direct imaging mode, which was used for the data pre sented here 



(all ta rgets in the range 6.5 < V < 13), achieves precision of ~0.1 millimags ( IRowe et al. 



2006b) 
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Photometric observations of AD Leo (V = 9.4) were obtained with MOST between 6-14 
March 2010, which corresponds to 3717-3726.5 in the MOST standard system of JD-J2000 
(JD - 2451545, the epoch of January 1, 2000). Our observations have a cadence of 60 seconds 
through a broad optical bandpass of ~3500-7000A. Guide stars were captured in 3-second 
exposures and stacked in groups of 10. 



Photometry was performed by the MOST pipeline (described in detail in iRowe et al. 



2006bl ). Several of the processing steps that are typically performed by the full pipeline, 
namely those that clean the light curve and attempt to remove spurious signals and outlier 
points, were not applied to the AD Leo data since they carried the risk of filtering out 
legitimate flare signals. The partially-processed light curve will hereafter be referred to as 
the "raw" light curve. 

The MOST observations, shown in the top panel of Figure [TJ produced a time series 
of 8,592 epochs over 9 days. The first 89 epochs, representing only 1% of the total, were 
followed by a 23- hour gap and were neglected from further analysis. Figured] (bottom) shows 
the raw time series after partial MOST pipeline processing and after a sinusoid was removed 
(described further in Section |2TT|) . 

Many exposures were missing from the raw light curve due to instrumental issues or 
suspected cosmic rays. These missing exposures are most noticeable near JD-JD2000 = 
3724.5, which is also shown in more detail in Figure [2] (top). The gaps in this period prevent 
us from successfully identifying and measuring all flares. We therefore additionally neglect 
the data in between JD-JD2000 = 3724.2364 and 3724.7151. Additional short gaps longer 
than three minutes are also removed from the overall monitoring time. Our final light curve 
contains 8383 points, corresponding to a total of 5.82 days of observations. 

The MOST satellite experiences phases of increased stray light either from the illumi- 



nated side of th e Earth flWa 



pass of its orbit (IRowe et al 



ker et al.ll2003l ) or the moon that occur regularly during each 
l2006bl ). These phases last for 43.2 minute intervals at a period 
of 101.4 minutes, and produce non-homogeneous intensity variations. Typically, this effect is 
removed from MOST light curves by folding the time series with the satellite's orbital period, 
calculating a running average, and correcting the entirety of the time series by this average. 
However, because the amplitude of the stray light modulation varies over time, there are 
over- or under-compensations in different parts of the time series, some of which can mimic 
a flare signal. Figure [2] (bottom) shows the most prominent of these fluctuations in the time 
series, between 2.3 and 4.6 days after the start of observations. The dimming ranges from 
~0.4 ADU pix -1 sec -1 below the median, to less than the rms scatter (<0.2 ADU pix -1 
sec -1 ). We therefore did not correct the lightcurve and instead account for the periodic 
dimming signature when searching for flares, as we discuss in detail in Section 12.21 
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2.1. Modeling and Removal of the Periodicity of the Light Curve 

The raw MOS T light curve of AD Leo was dominated by a large, regular modulation that 



we int erpret to be caused by a star spot rotating into and out of view, as in iPettersen et al. 



(119841 ) . The frequency of the modulation, co, was obtained using a Lomb-normalized pe- 
riodogram and is shown in Figure [3J We note that the orbital period of the spacecraft 
(cj = 14.2dy _1 ) does not produce a noticeable peak. We adopt the FWHM of a Gaussian 
fit to the periodogram peak as the uncertainty in the period. The frequency of 0.449±0.063 
dy _1 correspond s to a rotational pe riod of 2. 23^^27 dy, m good agreement with the 2.24 dy 



period found by iMorin et al.l ( 120081 ). 



Identifying flares in the light curve is best accomplished when any signatures not caused 
by flaring have been removed. An empirical model, F(t), of the form 

F(t) = - sm[2nu}(t - B)} + Ct + D (2) 

was fit to and then subtracted from the light curve. The parameters A through D (0.5161, 
0.9007, 0.0406, 0.2053) were determined by a least-squares fit of the model to the light curve, 
using the frequency (u) determined above. 

The linear term accounts for the small observed decrease in flux over the course of the 
observations that may be caused by a much longer time-scale modulation or by instrumental 
effects. In either case, it affects neither our period determination nor our flare identification 
and can therefore safely be removed. 

Figure [T] (bottom) shows the resulting flattened light curve with the empirical model 
removed. The small-scale structure remaining within the light curve is due to systematic 
noise (such as stray light from the illuminated side of the Earth - Section [2]) and intrinsic 
stellar variability, such as candidate flares. 



2.2. The Flare Identification Algorithm 



Fl are ca ndidates were identified using the Flare Identification Algorithm (FIA) of [Hilton 



(120111 . 120121 ). Briefly, at each epoch, i, the running local median and standard deviation 
(viocai) are computed using the 31 epochs nearest in time to % that have not been flagged as 
candidate flares. The procedure consists of iteratively identifying consecutive epochs that 
fall sufficiently far above the local median, flagging them as flare candidates, masking them 
so that they are not included in the determination of the local standard deviation, and 
repeating those steps on the masked light curve until new candidates are no longer found. 
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The criteria that were used to identify flare candidates required > 3 epochs above 2.5ai oca i, 
at least one of which is above 3.5<7z oca /. Flare candidate start and stop times were determined 
by identifying the last epoch prior to - and the first epoch after- the flare peak in which the 
running mean of 5 epochs was less than 0.5cr toca j. 

As discussed in £j2j some portions of the light curve contained periodic decreases in flux 
for which corrections were not applied. The FIA computes a running standard deviation that 
would become larger in the sections with decreased flux, thereby increasing the minimum 
flux threshold that the light curve would have to reach to be considered a flare. Because 
the increase in ai oca i is caused by systematic effects and not a decrease in the photometric 
precision, during the periods when the dimming occurs, we mask the dimmed regions in our 
computation of the running local median and standard deviation, in effect treating them the 
same as the flare epochs that are also removed. Although this masking allows us to detect 
flares near these regions, flares that occur during or near these periods must be slightly larger 
than other, smaller flares in order to be detected. For the light curve, the difference between 
the median values during the orbital dimming periods and outside of the orbital periods is 
less than 0.01 ADU pix -1 sec -1 . For the period shown in the bottom panel of Figure [2J the 
total time when the light curve is more than 0.1 ADU pix -1 sec -1 below the median is 2.6 
hrs. This is less than 2% of the total monitoring time and has negligible effect on the final 
FFD. 

Once the FIA was applied, the results were visually inspected and the flare candidates 
that did not resemble the classical flare shape of a sharp rise with an exponential decay were 
rejected. The classic flare is shown in Figure HI In one case, shown in Figure [2] (top), a 
candidate flare with a sharp rise and a portion of the initial decay phase occurred at the 
beginning of a period that was neglected because of missing data. We do not include this 
flare in the majority of our analysis, although we we will mention it briefly when we consider 
the flare wait time distribution (see Section I3~T1) . 

3. Results 

When applied to the MOST time series, the FIA resulted in 24 flare candidates, of which 
19 passed visual inspection. These flare candidates are shown in shaded regions in Figure 
while the detected but rejected flares are shown as hatched regions. One of the rejected 
flares is the truncated flare shown in Figure [2] (top) that occurred during a neglected portion 
of the light curve. The other rejected candidates occurred during the continuous portion of 
the light curve, but did not show the typical profile of stellar flares. They may comprise a 
series of short, complex flares that were not well-resolved, or they may be epochs that are 
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compromised by data loss, poor photometry, cosmic rays, etc. 

Of the 19 confirmed flares, the largest (Figure H]) reaches its maximum at 6 ADU pixr 1 
sec -1 above the median and lasts for 3.97 hrs. The peak was over 200 times the <Ji OC ai value, 
with a total energy release of logE = 33.5 ergs. The smallest flare was roughly 0.4 ADU 
pix -1 sec -1 above the mean, lasting for 0.12 hrs (~7 minutes). The peak was 10 times the 
ciocai value, with a total energy release of logE = 31.6 ergs. 

To compute the flare energy, we first convert the instrumental ADUs to phys ical units 



throug h spectrophotometry. We use a flux-calibrated spectrum of AD Leo from iKowalski 



( 120121 ) obtained with the Dual-Imaging Spectrograph on the ARC 3.5m telescope at Apache 
Point Observatory, and convolve it with the MOST bandpass. We then calculate the dimen- 
sionless flux from the flare, normalized by the quiescent flux of the star, as 

p /^s _ F Flare (t) + F _ ^ _ F Flare (t) 

Fq Fq 



wher e P is the quiescent flux. We then compute the equivalent duration, P (IGershberg 



19721 ). which is defined as the amount of time that it would take the star, in its quiescent 
state, to release the same amount of energy released during the flare and is simply the time 
integral of Ff. 

P= [F f (t)dt (4) 



The flare energy is the quiescent luminosity of the star, q, multiplied by the equivalent 
duration of the flare, P: 

E flare = Q X P (5) 



We stress that these values are filter dependent, such that flare energies measured in different 
filters are not necessarily equal. Finally, the cumulative FFD is found from the ordered 
list of flare equivalent durations, or Pj for the ith flare. The frequency of a flare with 
P > Pi is % (which is the number of flares larger than Pj) divided by the effective monitoring 
time. Because the MOST photometry has similar precision throughout the light curve, i is 
essentially divided by the cumulative amount of time where there are observations. Note that 
the total monitoring time does account for gaps in the otherwise continuous observations. 

We find that qMOST = 6.96 x 10 30 erg s" 1 using the flux determinations described above 
and a distance of 4.9 pc for AD Leo (jPerryman fc ESAlll997l ). However, because the MOST 
photometry is dominated by red light, and previous flare studies have used a narrower and 
bluer filter, such as Johnson U, the flare energies cannot be directly compared. Therefore we 
compare the frequencies, normalizing the equivalent durations to the largest value of log P, 
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allowing us to compare the values of f3 from Equation HJ The flare energies are used for the 
analysis of flare timing (see Section l3TTj) . 

The normalized FFD from our observations is sh own in Figure E] ( black line and points), 
along with the n ormalized FFDs for AD Leo from iPettersen et al.l (11984 . 115 flares) and 
Lacy et al.l (119761 . 9 flares). It should be noted that the fits are weighted towards high-energy 
flares, since all of these studies are affected by detec tion thresholds at the low energy end. We 
find = — 0.68 ± 0.16 , which is in agreement with lLacy et al.l (119761 . /3 = —0.82 ± 0.27) and 
Pettersen et al.l (11984 . f3 = —0.62 ± 0.09). We compute uncertainties in the same manner as 
these previous studies, with erg = where n is the number of flares. The general agreement 
is remarkable, especially considering that the other studies were performed in the Johnson 



U filter. Flares relea se more energy at bluer wavelengths (e.g., lHawley &: Pettersen! 11991 



Kowalski et al.l 120101 ). and AD Leo has much less quiescent flux in U than in the broader 
MOST filter. The broad filter means that despite the sensitivity of MOST, only relatively 
large flares are detected, and the overall flare rate is lower than it is in the U filter. 



3.1. Flare Timing 

The MOST light curve is nearly continuous over several days, with enough precision to 
detect photometric variations on timescales of a day or longer. We are thus able to investigate 
when flares occur relative to other flares and relative to the (presumed) rotational phase. 
These studies cannot be easily done from the ground due to diurnal observing interruptions. 

We wish to use the timing and energies of flares as probes of the distribution of magneti- 
cally active regions on the stellar surface. In this paper, we adopt a simple model that can be 
tested with the observations. We assume that all of the magnetic energy that is generated by 
the stellar dynamo is stored in active regions during quiescence, and is only released during 
flares. We also assume a single, large active region which stores the magnetic energy. This 
assumption is supported by the strong periodic signature in the MOST light curve, which 
we interpret as the rotational modulati on of a large active region. It is further supported by 



the spectropolarimetric observations of iMorin et al.l ( 120081 ). who report that one large polar 



region, viewed nearly pole-on, is the dominant magnetic feature on AD Leo. In our model, 
then, the wait time between a flare and the proceeding flare should be proportional to the 
flare energy. Furthermore, flares should be more likely to be observed when the active region 
is visible, which corresponds to the troughs of the light curve assuming that the active region 
is dark in the MOST filter. 



We plot the flare energies against the time since the previous flare in the top panel of 
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Figure [7J The bottom panel displays the number of flares in each 1 hour bin of wait time. 
The longest wait time between flares does occur prior to the largest flare in the set. The 
next longest delay is 18.95 hrs, leading to the flare candidate beginning at 3724.86 dy. While 
this flare candidate is not notably large, it must be noted that this is the first flare to occur 
after the neglected portion of the light curve. At least one flare appears to have occurred 
during this period (the truncated flare shown in Figure [21 top). If we consider this to be 
the most recent flare, the delay of 18.95 hrs is closer to 13.59 hrs, which is marked by the 
dashed line. In any case, the data in Figure [7J do not show a strong correlation between wait 
time and flare energy. In fact, some of the largest flares in the sample occur after very short 
delay times. 

In the top panel of Figure El we plot the flare energies against rotational phase, where a 
phase of 0.5 is the light curve minima. The bottom panel shows the number of flares in each 
rotational phase bin of size 0.1. The flares do not occur preferentially closer to minimum 
light, nor do the flares that occur near minimum light have larger energies. This indicates 
that either 1) any difference between flare occurrence at minimum light and maximum light 
is not detectable in our small sample, implying that the active region is large (or nearly pole 
on), such that even away from minimum light, flares that occur from the active region are 
visible; or 2) there are multiple active regions at a variety of longitudes that host flares, and 
while the rotational modulation is caused by the largest active region, flares may come from 
any of them. 

Taken together, the distribution of flare wait times and the phase distribution of flare 
occurrence do not support a non-circumpolar single spot model, however, this result may 
also be due to low number statistics. Future studies involving a much larger sample of flares 
from continuous monitoring, such as those that can be obtained with Kepler short-cadence 
(~ 1 minute) data, will provide more robust statistics of these two measurements that may 
allow us to reconsider this model, as well as consider more complex models. 



4. Conclusion 



We use over 8,500 individual MOST observation s of AD Leo t o find a rotational period 



of 2.23 q 27 dy. Using the automated, iterative FIA of lHiltonl (120121 ). we identify 24 flare can- 



didates, of which 19 pass visual inspection. The FFD derived from these flares is consistent 
with previous studies, although it suffers from low number statistics and the observations 
were obtained with a much broader and redder filter. 



The continuous monitoring provided by MOST allowed for a previously-unattainable 
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wait-time analysis. Although the long waiting period leading up to the most energetic flare 
in this study suggests the possibility of a relationship between flare delay times and the 
occurrence of large, high-energy flares, no meaningful relationship is apparent. Continuous 
monitoring and high precision photometry also allowed for the investigation of flare occur- 
rence with rotational phase. We find no evidence of a relationship, although low-number 
statistics prevent us from ruling out the existence of such a relationship. 

We compare these observational results to those expected from a simple model of mag- 
netic energy storage and release. In the model, magnetic energy is produced at a constant 
rate and stored in a single large magnetic spot, which acts as a reservoir. When the reservoir 
fills, magnetic energy is released in the form of flares. The model predicts there to be a rela- 
tionship between the wait-time of flares and flare energy, as well as between flare occurrence 
and rotational phase. The model predictions are not supported by the observations. This 
discrepancy could be explained if there are multiple reservoirs of magnetic energy or if the 
reservoir is not req uired to be full for a flare to occur. One possible scenario, supported by 



Morin et al.l (120081 ). is of a large magnetic spot viewed nearly pole on. This scenario predicts 
there to be a weak or nonexistent relationship between rotational phase and flare occurrence, 
but still predicts there to be a relationship between flare energy and wait-time. 

Further studies with longer monitoring times and enhanced photometric precision are 
necessary. The Kepler mission, in its short-cadence mode, allows for additional studies, the 
results of which will be presented in future papers. 
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Fig. 1. — Top: The raw MOST light curve of AD Leo. The gray sinusoid is the calculated 
rotational period. Bottom: Residual model-subtracted light curve with stray light modula- 
tion and flare candidates. Epochs prior to 3718.27 JD-JD2000 dy have been removed due to 
insufficient observations. 
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Fig. 2. — Top: An enlarged section of the light curve, showing the truncated flare between 
3724.2 and 3724.4 dy (hatched region). This flare is not counted amongst the final flare 
candidates due to insufficient data. Bottom: Modulation of the light curve due to stray 
light, from 2.3 - 4.6 days after 3717.1084 JD-JD2000 dy. The period of this modulation 
corresponds to the 101.4 minute orbital period of MOST. Epochs of dimming due to this 
modulation were masked when calculating the FIA. 
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Fig. 3. — Results of the Lomb-normalized periodogram showing the power corresponding 
to the peak power of the periodogram at a rotational frequency of u = 0.449 ± 0.063 dy -1 
(dashed line ) . Th is period agrees well with the 2.2399 dy {u = 0.4464 dy -1 ) period of 
Morin et al.l (120081 ) . The orbital period of the spacecraft (u) = 14.2dy— 1) does not produce 
a significant peak. 
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Fig. 4. — An enlarged section of the light curve showing the largest flare candidate in the time 
series. This flare lasts for 3.97 hours, beginning at 3721.97 dy and returning to quiescence 
at 3722.05 dy. The peak of this flare is at 28.83 ADU pix -1 sec -1 representing an increase 
of 28% over the average quiescent value of the time series. 
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Fig. 5. — The model-subtracted light curve with stray light modulation. Solid shaded regions 
identify the 19 flares that passed visual inspection out of the original 24 candidates. The 
rejected flare candidates are shown as hatched regions. 



-18 - 




Fig. 6. — Cumulative flare frequency distributions (FFD), normalized to the maximum 
equivalent durations (P) for each distribution. The FFDs have been normalized to facilitate 
comparison between the observations take in the MOST filter and those taken in the U filter. 
The MOST obs ervations (solid l ine) sh ow less frequent flaring, but have a similar j3 value as 
the FFDs from iPettersen et all (11984 dotted line) and lLacy et al.l (119761 dashed line) . 
The uncertainty in j3 for the MOST observations is shown as the shaded region. 
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Fig. 7. — Top: Flare energy is not strongly correlated with the time since the previous flare. 
Although the largest flare in the sample also has the longest wait time, there are several large 
flares with very short wait times. The second longest wait time occurred after the truncated 
and rejected flare in Figure [2] (top). Including the rejected flare reduces the wait time to 
~12 hours, as shown by the dashed line. Bottom: The distribution of wait times between 
flares in bins of 1 hour, shows that most flares occur within a few hours of the previous flare. 
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Fig. 8. — Top: Flare energy E vs. rotational phase (f) rot , where 0.5 is minimum light. 
Bottom: Distribution of the number of flares in phase bins of 0.1. We might expect to see 
more flares during minimum light, when the active region causing the modulation is most 
visible. However, this is not apparent in our current observations. 



